Bottomonia in the Quark-Gluon Plasma and their Production at RHIC and LHC 
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We study the production of bottomonium states in heavy-ion reactions at collider energies avail- 
able at RHIC and LHC. We employ an earlier constructed rate equation approach which accounts 
for both suppression and regeneration mechanisms in the quark-gluon plasma (QGP) and hadroniza- 
tion phases of the evolving thermal medium. Our previous predictions utilizing two limiting cases 
of strong and weak bottomonium binding in the QGP are updated by (i) checking the compatibility 
of the pertinent spectral functions with lattice-QCD results for euclidean correlators, (ii) adapting 
the initial conditions of the rate equation by updating bottom-related input cross sections and the 
charged-particle multiplicity of the fireball, and (iii) converting our calculations into observables 
as recently measured by the STAR and CMS experiments. Our main findings are a preference for 
strong T binding as well as a significant regeneration component at the LHC. 



I. INTRODUCTION 

Quarkonia by now have a 25-year history of being uti- 
lized as a probe of the Quark-Gluon Plasma (QGP) in ul- 
trarelativistic heavy- ion collisions (URHICs), cf. Refs. p]- 
Q for recent reviews. The majority of the studies to date 
have been devoted to charmonium properties and their 
observable consequences in experiment, with few excep- 
tions for bottomonia Recently, experiments at the 
Large Hadron Collider (LHC) [10L lllj] and the Relativistic 
Heavy-Ion Collider (RHIC) [12| have measured T mesons 
in URHICs for the first time, and subsequent phenomeno- 
logical analyses have commenced [13, [ljj], focusing on 
suppression mechanisms in Pb-Pb( v / s=2.76 ATcV) in 
comparison to CMS data. For charmonia, regeneration 
effects due to coalescence of c and c quarks have compli- 
cated the interpretation of observables appreciably. Since 
the cross section for open bottom is much smaller than 
for open charm, regeneration of bottomonia has been ar- 
gued to be negligible. This would render bottomonium 
production in URHICs a more pristine probe of quarko- 
nium dissolution in the QGP. However, quantitative esti- 
mates @, Q indicate that this may not be the case. The 
magnitude of the regeneration contribution depends on 
a subtle interplay of the masses of the open- and hidden- 
bottom states in a system with fixed bb content. In ad- 
dition, the ratio of hidden- to open-bottom states in el- 
ementary hadronic collisions is typically of order 0.1% - 
and thus much smaller than in the charm sector where it 
is around 1% -, implying that even small contributions to 
bottomonium regeneration can be significant relative to 
primordial production. Furthermore, a consistent inter- 
pretation of bottomonium production in URHICs should 
encompass both LHC and RHIC data. 

The above aspects were addressed in previous work Q 



where predictions for the nuclear modification factor, 
Raa , of inclusive T production at RHIC and LHC have 
been presented. Unfortunately, these calculations could 
not be directly mapped to the now available data: at 
RHIC the Raa for the sum of IS, 2S and 3S states has 
been measured, and at LHC the current data are taken at 
half of the design energy, -^5=2. 76 ATeV. Furthermore, 
the availability of bulk-hadron observables in Pb-Pb col- 
lisions at the LHC helps to constrain the fireball model 
in terms of its total entropy and transverse expansion 
which govern the temperature evolution, thus affecting 
bottomonium suppression and regeneration. Calcula- 
tions of the latter also benefit from an updated input 
for the bottom(onium) cross sections (if only suppres- 
sion is considered, this input cancels out in the suppres- 
sion factor). We also improve the theoretical input to the 
earlier calculations by applying theoretical constraints on 
the in-medium bottomonium properties (binding energy 
and dissociation width) in terms of euclidean correlators 
computed in recent years within thermal lattice-QCD 
(1QCD) [H-[l3l- InRef. 0, two scenarios for in-medium 
bottomonia were discussed, i.e., one with vacuum and 
one with in-medium binding energies estimated from a 
screened Cornell potential [18| . With the now available 
lattice data, we can refine the pertinent spectral func- 
tions following the approach outlined in Rcf . [lj| . 

Our paper is organized as follows. In Sec. [H] we re- 
visit the spectral properties of bottomonia in an equili- 
brated QGP within two limiting scenarios, "strong" and 
"weak" binding, and check the resulting euclidean cor- 
relators against 1QCD data. In Sec. IIIII we briefly recall 
the kinetic rate-equation framework in a thermal medium 
and discuss its (suitably updated) input quantities. In 
Sec. IIVI we present the updated results on bottomoniun 
production at RHIC and LHC, focusing on its central- 
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ity dependence in comparison to STAR and CMS data. 
We evaluate our results and put them into context with 
existing calculations. In Sec. [V] we summarize and con- 
clude. 



II. BOTTOMONIUM SPECTRAL FUNCTIONS 
IN THE QUARK-GLUON PLASMA 

Following Ref. Q we attempt to bracket the properties 
of bottomonia in the QGP by defining two scenarios for 
the temperature dependence of their binding energy, 

e B (T) = 2m b (T)-m Y (1) 

(m&(T): in- medium bottom mass; my- bottomonium 
bound-state mass with Y = T(1S), T(2S), X&(1P), . . . ). 
These are: 

1. Weak Binding Scenario (WBS): bottomonium 
bound-state energies are significantly reduced in 
the QGP, relative to the vacuum state, adopting 
the screened Cornell-potential results of Ref. fl8j |. 
The bound state masses are kept at their vacuum 
value which is motivated by the approximate con- 
stancy of their in-medium correlator ratios as found 
in 1QCD (further details below). This, in turn, re- 
quires the in-medium bottom quark-threshold to 
decrease with increasing temperature (which is not 
inconsistent with independent 1QCD computations 
of the heavy-quark free energy) . 

2. Strong Binding Scenario (SBS): the temperature 
dependence is removed from Eq. (TTJ), implying 
bound states to remain at their vacuum mass and 
binding energy throughout the QGP, with a con- 
stant open-bottom threshold at 2mf, = 2ms = 
10.56 GeV. 

Figure [1] shows the temperature dependences of cb for 
the ground state T, for T' and x&(lP), within these 
two scenarios. The constant values in the SBS are to 
be contrasted with factor of at least 2 reduction at 
T c =180MeV (and rapidly decreasing further with T) 
in the WBS. Strong- and weak-binding scenarios also 
emerge in potential-model calculations, due to current 
uncertainties in defining the in-medium interaction ker- 
nel, usually bracketed by the internal and free energies 
of the QQ, see, e.g., Ref. (20[. In the present context we 
arc not (yet) aiming at a quantitative implementation of 
potential-model results, but rather at exploring the sen- 
sitivity of heavy-ion data to Y spectral properties within 
the two cases defined above. 

The next step is the evaluation of the Y dissociation 
widths. For a Coulombic bound state with large binding 
energy, es 3> Aqcd, the leading break-up mechanism 
has been identified and calculated more than 30 years 
ago as the gluo-dissociation process, g + Q —> Q + Q 
(Q: quarkonium) plj . The corresponding cross section is 
peaked around incoming gluon energies of fcp = (10/7)es. 
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FIG. 1: Temperature dependence of the binding energies, 
£b(T), for T (solid lines), T' (dotted lines) and Xb (dashed 
lines) in the Strong- and Weak-Binding scenarios (upper and 
lower 3 lines, respectively). 



For weakly bound quarkonia, < T, this mechanism 
becomes inefficient and is supplanted by the quasifrec 
break-up, p + Q—>p + Q + Q(p = g, q, q) |22{ . This 
picture has been well established by now. For exam- 
ple, in cffectivc-ficld-thcory approaches to quarkonia, the 
gluo-dissociation corresponds to the so-called "singlet-to- 
octet" transitions [2^, [24| while the quasifree mechanism 
corresponds to Landau damping in the gluon-cxchangc 
potential [23|,[2i| . Accordingly, we will treat the SBS with 
the gluo-dissociation and the WBS with the quasifree 
mechanism. The pertinent dissociation rates (inelastic 
widths) can be expressed in terms of an underlying cross 
section, a Y lss , convoluted over the thermal parton distri- 
bution functions, f p , as 

T*T{q;T) EE Ty'^T) 

= E / -^f P (uk-,T)a Y iss ( S )v rci (2) 
p=q,q,g 

where u)k denotes the on-shell energy of a parton with 
3-momentum k, and s = (k^ +q^) 2 is the total center- 
of-mass energy squared in the collision with a Y of 3- 
momcntum q (v Te i is their relative velocity). The bot- 
tomonium lifetimes arc plotted in Fig. [2] for q = 0, for 
T, T' and Xb i n the weak- and strong-binding scenar- 
ios (left and right panel, respectively). As in Ref. Q 
quasifrec dissociation is calculated with a s ~0.26, while 
for gluo-dissociation a Coulomb-like binding for T(1S) 
with es=l.l GeV translates into an effective coupling of 
a s ~0.65. Within the WBS, quasifrec dissociation in- 
duces little difference between the T' and \b widths, due 
to their almost identical (small) binding energies. The 
T(1S) lifetime is about an order of magnitude longer than 
for the excited states at T c , but quickly approaches their 
lifetimes since eg decreases to small values at higher tem- 
perature. In the SBS gluo-dissociation produces a larger 
spread in the lifetimes of ground and excited states. At 
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FIG. 2: Bottomonium lifetimes in the QGP for the two binding scenarios defined in the text; left panel: WBS with quasifree 
dissociation; right: SBS with gluo-dissociation; solid lines: T, dashed lines: T', dotted lines: \b- 



T c , the T(1S) lifetime is by almost two orders of magni- 
tude larger than its WBS counterpart, but it comes close 
to the latter at high temperatures, T ~ 1 GeV. However, 
it stays well above the lifetimes of the excited states even 
at high T. Close to T c . the lifetimes of the excited states 
in the SBS are quite comparable to their WBS coun- 
terparts, even dropping below the latter for T>250MeV 
or so. The reason for that is the larger coupling in the 
SBS together with a favorable overlap of the typical ther- 
mal gluon energies, uik — 3T, with the maximum in the 
gluo-dissociation cross sections around Wfc^es~600 McV. 
Thus, for the phenomenologically relevant temperatures, 
especially at RHIC, one expects the suppression factors 
of the excited states to be similar in the SBS and WBS. 

With the masses and widths (real and imaginary parts 
of the inverse propagator) of the Y bound states de- 
termined we proceed to compute the pertinent spectral 
functions. We follow Ref. [19| for the charmonium case 
by adopting an ansatz consisting of a relativistic Breit- 
Wigner (RBW) for the bound-state part and a nonper- 
turbative continuum for energies above the bb threshold. 
We neglect the excited states of Y. One has 



ov(w) = A-] 



2ui 



ZyuiTyiT) 



tt (u; 2 -m Y ) 2 +w 2 IY(T) 2 



(3) 



where -ZV C =3 is the number of colors, Ay is given by the 
wave function overlap at the origin in vacuum, Zy char- 
acterizes its modification in medium ("polcstrength"), 
s(T) = 2mb(T) marks the in-medium bb threshold follow- 
ing Eq. (JXJ) , and the coefficients (a, b) arise from the Dirac 
structure specific to the hadronic channel under consid- 
eration, e.g., (2,1) for the vector channel. The prefactor 
By = 2 augments the perturbative expression of the con- 
tinuum to account for bb rescattering effects and is esti- 
mated from the T- matrix calculations of Ref. [2(|. We 
approximate the total width for the bound state with its 



dissociation rate (inelastic width). Figure [3] shows the in- 
medium T spectral function for two temperatures in the 
WBS and SBS. One clearly recognizes the much larger 
width and reduced threshold in the WBS compared to 
the SBS. In the latter the bb threshold remains constant 
at ^fs = 10.56 GeV. 

In a final step we check the above spectral function 
against results from thermal 1QCD where the bottomo- 
nium correlators, G(r), are computed as a function of 
cuclidcan time r. The euclidean correlator is related to 
the spectral function as 



G Q (r;T) = / dcu a a (cj;T) K(u,t;T) , 



(4) 



where a specifies the hadronic channel, and the kernel 
reads 



K(u,r,T) 



cosh [lo(t- 1/2T)] 
sinh[w/2T] 



(5) 



To better exhibit the medium effects, it is common to 
consider the correlator ratio, 



Rg(t;T) 



G a (r, T) 
G r a ec (r,T) 



(6) 



where the denominator is evaluated with the same ker- 
nel but with a spectral function evaluated at a low (or 
zero) temperature, T*, where no significant medium ef- 
fects are expected. To compute the euclidean correlator 
ratios for our in-medium spectral functions we define a 
vacuum spectral function using Eq. (J3J) in the narrow- 
width limit and with Zy=l for the RBW part and with 
^/s = 2iiiB = 10.56 GeV for the continuum (as for the 
in-medium spectral function, we only include the ground 
state). Our results for Rq(t,T) arc displayed in Fig. 2] 
In the strong-binding scenario (featuring vacuum bind- 
ing energies and bb threshold) , the only medium effect is 
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FIG. 3: Spectral functions of the T calculated from Eq. Q in a QGP of temperatures 2T C (solid lines) and 3T C (dashed lines); 
left panel: weak-binding scenario, right panel: strong binding scenario. 
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FIG. 4: Euclidean correlator ratios following from the in-medium T spectral functions shown in Fig. [3] at various temperatures 
for the weak-binding (left panel) and strong-binding scenario (right panel). 



the increasing decay width, which, however, has very lit- 
tle impact on the correlator ratio so that the latter stays 
within 2% of one up to the highest temperature consid- 
ered (4T C ), cf. right panel of Fig. @] In the weak-binding 
scenario, the bound-state mass is also constant, but the 
threshold energy drops appreciably, which by itself would 
imply an increase of R(t) above one. However, a weaken- 
ing bound state is not only characterized by a reduction 
in eg but also by a reduction of its polestrength (see, 
e.g., Ref. which is represented by the Zy-factor in 

the RBW. We have tuned this quantity to recover R(t) 
at around one resulting in the curves shown in the left 
panel of Fig. [4] The deviations from one are larger in 
the WBS compared to the SBS, comparable to what has 
been found for charmonia with the same ansatz [l9j . The 
SBS, therefore, appears to be a more natural realization 
of the near-constancy of R(t) ~ 1 at the few-percent 
level as found in thermal 1QCD [l5l - [l7j . This insight is 
one of the main motivations for re-evaluating the previ- 



ous phenomcnological analysis of Ref. [7[ where the in- 
medium binding scenario with quasifree dissociation was 
believed to be more realistic for bottomonia in the QGP. 
To our knowledge, this has not been done in any other 
phenomenological application thus far. 



III. KINETIC-THEORY MODEL 

In this section we first set up the rate-equation frame- 
work for the Y states in the expanding QGP (Sec. IHI Al) . 
followed by specifying the inputs and initial conditions 
for their numerical evaluation, specifically 

(i) the temperature and volume evolution of the ther- 
mal bulk medium (Sec. IIIIBj) . and 

(ii) the initial conditions for the various Y states and 
the number of bb pairs in the system (Sec. IIIICp . 
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A. Rate Equation 

To simulate the evolution of the bottomonium abun- 
dances in URHICs at RHIC and LHC we adopt the rate- 
equation approach of Ref. (2(| (also used in Ref. Q). It 
describes the approach of the number, Ny, of bottomo- 
nium state Y, toward equilibrium as 



d/Vy 
~d^ 



-pdiss 
_i Y 



(T) [Ny 



Ny q (T)} 



(7) 



The key quantities in this equation are the inelastic re- 
action rate, Fy ISS , and the equilibrium limit, N Y q , which 
directly relate to the spectral properties discussed in the 
previous section. The inelastic reaction rate would be 
all that is needed if only suppression effects were consid- 
ered, governed by the first (loss) term on the right-hand- 
side of Eq. 0. Regeneration processes are accounted 
though the gain term, r Y iss (T)N Y q (T). Its form is dic- 
tated by the detailed balance between the dissociation 
and regeneration processes. Apart from the inelastic re- 
action rate, the gain term is controlled by the equilibrium 
limit which is computed as follows. For a given collision 
energy, system and ccntrality, the total number of bb pairs 
is assumed to be constant throughout the evolution of 
the expanding fireball. Assuming relative chemical equi- 
librium between all available states containing bottom 
quarks at given temperature and volume of the system, 
the bb number is matched to the equilibrium numbers of 
bottom states by the condition 



N, 



bb 



1 h(N op ) 
2 N ° P MN^) +Nhid 



(8) 



by means of a 6-quark fugacity factor, 7& = 75. The equi- 
librium open- and hidden-bottom numbers in the QGP 
are given by 



7b Vfb d b 



d 3 P i 
(2tt) 3 



(9) 



hid 



V FB J2 d Y I 



f Y (p-:T), (10) 
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FIG. 5: Equilibrium numbers of T vs. temperature at RHIC 
(^/snn = 0.2 TeV, lower two lines) and LHC (^snn = 
2.76 TeV, upper two lines) using equation (|10[) with bottom- 
production cross sections as specified in Sec. lIHCl below, The 
solid lines correspond to the weak-binding scenario (with in 
medium 6-quark mass) and the dashed lines to the strong- 
binding scenario (with vacuum &-quark mass). 



In heavy-ion collisions, however, 6-quarks are not ex- 
pected to kinetically equilibrate |27| . Bottom-quark spec- 
tra which arc harder than in equilibrium imply less phase- 
space overlap between bb pairs and thus a reduced Y re- 
generation. To account for this we follow Ref. [28| by 
introducing a relaxation-time factor 



TZ(t) = 1 - exp 1 



dr' 



' eq 



(ii) 



(2*) 



whose timcscale is set by the kinetic relaxation time, r eq , 
of b quarks in the QGP. In our calculations it is taken 
as T oq ~ lOfm/c around T ~ 2T C , in a range consis- 
tent with microscopic T-matrix calculations of bottom 
diffusion (2(| (with these diffusion coefficients a fair de- 
scription of the observed open heavy-flavor suppression 
and elliptic flow at RHIC is possible [29(). 



and Vfb is the fireball volume at given temperature, 
T(r). The equilibrium number densities are simply those 
of b quarks (with spin-color and particle-antiparticle de- 
generacy db=6x2) and bottomonium states (summed 
over including their spin degeneracies, dy). They are cal- 
culated using the same in-medium masses for open- and 
hidden-bottom states as in the spectral function analysis 
in Sec. |TT] for both the SBS and WBS, respectively. At 
fixed N b i a large 6-quark mass leads to a large jb which 
in turn results in a large equilibrium number, Ny. This 
is reflected in the temperature dependence of the T equi- 
librium numbers which are much larger in the strong- 
than in the weak-binding scenario, cf. Fig. [5] The ratio 
of modified Bcssel functions, Ii(N op )/Io(N op ), in Eq. (|8J) 
enforces exact conservation of the net b number in each 
event (canonical ensemble). 



B. Fireball Model 

Our fireball model has been updated compared to 
Ref. Q by matching it to our recent charmonium applica- 
tions at RHIC [TJ| and LHC [3(| . We assume an isentrop- 
ically expanding isotropic firccylindcr whose eigcnvolume 
is approximated by the ansatz 



Vfb(t) = (z + v z t) tt [ r + ^a±r' 



(12) 



with initial transverse radius, tq, longitudinal velocity 
v z =lAc and transverse acceleration aj_=0.1c 2 /fm. The 
latter is larger than in Ref. Q to better reproduce mea- 
sured bulk-hadron p t spectra and leads to a slight re- 
duction in the QGP and mixed-phase lifetime. At each 
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FIG. 6: Temperature evolution of the fireball model at RHIC 
(solid line) and LHC (dotted line). At RHIC the evolution 
is given for central Au-Au collisions at ^/sjvjv=0.2 TeV and 
a QGP formation time to = 0.6 fm/c; at LHC the evolution 
is for central Pb-Pb collisions at y / ijvjv=2.76 TeV and to = 
0.20 fm/c. 



centrality and collision energy of a heavy-ion collision, 
the total entropy, S, is chosen to reproduce the ob- 
served charged-hadron multiplicity (computed with a 
hadron resonance gas at a chemical freezeout temper- 
ature, T ch =T c =180MeV). With S = 10000(22000) for 
iVpart=375 nucleon participants at RHIC (LHC) one ob- 
tains ciJV c h/(i2/=800(1750). The entropy density, s(r) = 
S/Vfb( t ), is used to determine the temperature, T(r), 
by means of a QGP quasiparticle equation of state for 
s(T). At T=T C a standard mixed-phase construction 
is employed. We have verified that the impact of the 
subsequent hadronic evolution (also in the mixed phase) 
is negligible for all Y states. The initial longitudi- 
nal length, zq = ToAy, controls the initial temperature 
(Ay = 1.8 is the rapidity coverage of the fireball). With 
the QGP formation time To = 0.6(0.2) fm/c one obtains 
T = 330(610) MeV at RHIC (LHC). The temperature 
evolution of the fireball for central collisions is displayed 
in FigjBJ with QGP and mixed phases indicated. 



C. Primordial Y and Bottom Yields 

The rate equation requires an initial condition for 
the number of bottomonia at the QGP formation time, 
Ny{ t o)- We determine these by using their produc- 
tion cross section in pp collisions and binary collision 
(N co n(b)) scaling in AA collisions. We use inelastic NN 
cross sections of a l ^=A2mb and 62 mb at y / s=0.2 and 
2.76 TeV, respectively. In addition, we allow for cold- 
nuclear-matter (CNM) suppression, collectively defined 
as being due to modifications of Ny prior to thermal- 
ization, most notably nuclear modifications of the par- 
ton distribution functions, absorption on passing-by pri- 
mordial nucleons (or even secondary particles), and the 



Quantity 


■^/s NN =U.'2 TeV 


Ref. 


y/s NN ='2.7b TeV 


Ref. 


[mb] 


42 


[32] 


62 


[32] 


ai^Ur [nb] 


6.6 


[33J* 


250 ± 30 


[34]* 


cr£+ r [nb] 


3.4 


[33J* 


128 


[35J* 


a-pp^r' [nb] 


2.2 


[33]* 


83 


[35J* 


a PP-^x b I nb ] 


7.1 


[33]* 


270 


[35J* 


[^b] 


3.2 ± 1.8 


[3JJ* 


142 ± 25 


[38]* 


Ry 


0.52 


[40] 


0.3 


[40] 


S(N P =375) 


10000 


[22] 


22000 


[30] 


t [fm/c] 


0.6 


[19] 


0.2 


[30] 



TABLE I: Summary of cross sections in pp collisions used 
to initialize the hidden- and open-bottom abundances of our 
rate equations in Au-Au at RHIC (second column) and Pb- 
Pb at LHC (fourth column). Uncertainties are quoted when 
quantitatively available; starred references require further ex- 
planation, given in the text. 



Cronin effect. For simplicity, we approximate the combi- 
nation of all CNM effects by a suppression factor, 



5 nuc = cxp [-p N <j ahs L(b)\ 



(13) 



with an effective nuclear absorption cross section, <7 a b s . 
The latter will be considered over a range of 1-3.1 mb 
at RHIC and 0-2 mb at LHC. The former is compatible 
with the STAR measurement of R dA u = 0.78±0.20±0.28 
for T(1S+2S+3S) at RHIC [HI, while CNM effects are 
expected to be reduced at LHC due to the much in- 
creased Lorentz contraction of the incoming nuclei. The 
other parameters in Eq. (|13p are the nuclear density, 
Pat=0.14 fm -3 , and the impact-parameter dependent 
path length, L(b), evaluated with a Glauber model for 
the nuclear overlap. 

Current experimental limitations prevent resolving 
of individual Y(nS) states at RHIC and are reported 
as a combined cross section for T(1S+2S-|-3S) in- 
cluding the T(nS) — > e + e~ branching ratios, i.e., 

J2n=i £>(nS)(r(nS) EH- The measured value turns out 
to be consistent with next-to-leading order (NLO) color- 
evaporation model (CEM) calculations (using MRST HO 
parton distribution functions) for the individual T pro- 
duction cross sections [33[ . We therefore adopt the latter 
in our calculations, which are close to the values used in 
previous work [3]. 

To estimate the bottomonium cross sections 
2.76 TeV we start from 
of T measured by 



at 



the inclu- 

CMS m, 



(employing 
Guided by 
to 2.76 TeV 



sive production 

a^ T (v^=7TeV,|y|<2)~(300±40)nb 
a fJ, + fi~ branching ratio of 
NLO calculations [42| we extrapolate 
by a factor of 0.5 and to full rapidity by a fac- 
tor of 1/0.6. For a rapidity window of Ay=1.8 
for one thermal fireball (^=0.3, see below), we 
obtain a^ T (v / s=2.76TeV,|j/|<0.9)~ (75±8)nb. 

Alternatively, we can estimate this quan- 
tity by taking the CMS midrapidity value of 
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aJ^ T (v/5=2J6TeV,|j/|<0.9)~(72±8)nb 
tent with the above estimate 



dapper /dy(yfs=7TeV) ~ (80±9)nb, downscale it 
by a factor of 2 and multiply by Ay=1.8 to obtain 

consis- 
These estimates 
are also compatible with the CDF value [43[ of 
dapper /dy(y/s=l. 8 TeV) ~ (27±1.5)nb, extrapolated to 
2.76 TeV with a factor of —4/3. Based on the inclusive 
T production, we fix the cross sections for direct T, T', 
and Xfc(lP) to maintain compatibility with the fecddown 
fractions given in Ref. [H| and listed in Tab. HU e.g., 
with a ca. 30% branching for T' — > TX one needs 
Cpp~yT' — g0pp->T to obtain a feeddown fraction of 
~10%, and likewise for the Xb states. 

The regeneration component requires the knowl- 
edge of the 66 cross section, which we also take 
from pp measurements together with binary-collision 
scaling. At RHIC the STAR measurement (36j 
gives cr*° t _ >bB ( % /i=0.2TeV) = (1.6±0.5)//b, where we 
combined statistical and systematic errors. This 
is slightly smaller than, but compatible with, the 
value of 2 fib used in Ref. Q. The recent PHENIX 
measurements finds a larger central value j37j . 
f7 t ° t _ >65 ( 1 /i=0.2TeV)=(3.2±1.8)/ib, which, however, 
is still compatible with the STAR datum within com- 
bined statistical and systematic errors [4f|. In our 
calculations below we will use a PHENIX value since it 
results in a T /bb ratio which is closer to the LHC value, 
but wc will also comment on using the STAR value. 
At LHC our baseline value is obtained from LHCb 



data at 7 TeV [34|, again downscalcd by a factor of 2. 



With R y =0.3 for one thermal fireball (see below), one 
then finds CT pp _ >fc6 (Vi=2.76TcV,|y|<0.9)~(43±8)^b. 
Alternatively, one can estimate this quantity by 
using the midrapidity values of (35±5) fib at 
7TeV [13 and (15±2) fib at 1.96 TeV [H to ex- 
trapolate to — (20±3)/ib at 2.76 TeV (per unit rapid- 
ity), which, upon multiplying with Ay=1.8, yields 
c r pp^ & b(%/s=2.76TcV,|y|<0.9) ~ (36±5) fib, consistent 
with the above value. Since nuclear shadowing is 
believed to be significant at LHC, and applicable to 
the initial bb production, a reduction factor of 0.75 
will be included in Pb-Pb collisions (note that for T 
states the shadowing effect is associated with the CNM 
suppression factor, S n uc, in Eq._ (fT3|) ). At RHIC, we 
assume no shadowing on initial 66 production. 

As indicated above, a thermal fireball of light particles 
covers about 1.8 units in rapidity; the total number of 
initially produced Y (or 66) states within the fireball are 
thus obtained as 



N x (t ) 



pp^fX 
incl 



Prompt T(ls) 


-51% 


T(ls) from X&(1P) decays 


-27% 


T(ls) from xt(2P) decays 


- 10% 


T(ls) from T(2S) decays 


~ 11% 


T(ls) from T(3S) decays 


- 1% 



TABLE II: The decomposition of the total (inclusive) T(ls) 
yield into the direct T and feeddown from excited states. The 
composition was determined from ^/sjvjv=39 GeV p-p colli- 
sions 0. 



magnitude of this contribution in pp collisions is given 
in Tab. UH reaching ca. 50% for the inclusive T yield. 
We assume this to carry over to the initial conditions 
for AA collisions and solve individual rate equations for 
direct T, T' and Xb] we do not distinguish individual Xb 
states so that their combined feeddown amounts to 37% 
in pp; the T' contributes 11%. The feeddown from the 
T" is considered negligible (1%). This decomposition of 
the total T yield is assumed to be independent of the 
collision energy and colliding systems. 



IV. COMPARISON TO DATA 

Utilizing the previously discussed binding scenarios 
and dissociation mechanisms (Sec. |TT|, implemented 
into the rate equation in a thermal-fireball background 
(Sees. IIII Al and IIII B[) . observables are calculated for 
comparison to recent data at RHIC and LHC in the form 
of the nuclear modification factor, 



RAA{N palt ) 



N AA 



(JVpart) 



AUi(JV, 



part ) 



(15) 



(T 



N coll (b) R y S n 



(14) 



The numerator, N AA , is the solution of the rate equa- 
tion (J7J) at chemical freezeout. Deviations from unity 
are induced by CNM effects, suppression and regenera- 
tion in heavy-ion collisions. All previous studies of R\a 
thus far show suppression to be the dominant effect; sig- 
nificance of regeneration in the total yield was only sug- 
gested in our previous calculation Q and in the statistical 
model Q where all primordial bottomonia are assumed 
to be dissociated. In the following figures, our calcu- 
lations for Raa are broken up into nuclear absorption, 
primordial, and regeneration components, a luxury not 
available to experiment, but instructive for interpreting 
the data. 



with X = T,T',x&,66; here, R y ~ 0.52(0.3) denotes 
the fraction of the total number of X contained in 
Ay=1.8 around midrapidity at RHIC (LHC). For X = bb, 
S' nuc =l(0.75) for RHIC (LHC). 

To properly account for the observed inclusive T yields, 
feeddown from higher bound states is computed. The 



A. RHIC 

As mentioned above our calculations at RHIC supple- 
ment previous work Q for the same collision system, 
with slightly updated fireball, spectral-function and in- 
put parameters. In Fig. [7] we summarize our results for 
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the i?AA(A r P art) in 0.2 ATcV Au-Au collisions, for the 
combined Y(1S+2S+3S) (including e + e _ branching ra- 
tios) [i^ compared to STAR data in the upper row, and 
for the direct T, T' and \b hi rows 2-4. The left panels, 
representing the weak-binding scenario, show substan- 
tial suppression, reaching a factor of ^2 for the ground 
state and more than 5 for the excited states in central 
Au-Au. Regeneration does not play any role due to a 
small equilibrium limit, recall Fig. [5] The combined 
i? J 4 J 4(lS+2S+3S) is in the lower parts of, but not in- 
consistent with, the STAR data. In the strong-binding 
scenario the suppression of the ground state is notice- 
ably less pronounced than in the WBS, being almost 
entirely due to CNM effects. On the other hand, it 
might seem surprising that the suppression of the ex- 
cited states is rather comparable in both scenarios. As 
discussed in Sec. UH this is due to the similar dissociation 
rates of Y(2S) and X&(1P) in SBS and WBS at moder- 
ate QGP temperatures (recall Fig. [5]), originating from 
the larger coupling and favorable dissociation kinemat- 
ics in the gluo-dissociation processes (thus compensating 
the larger binding energy). This, in particular, puts into 
question a "sequential suppression pattern" (or "ther- 
mometer") which is often discussed without taking into 
account the finite width of the bound states. The weaker 
suppression in the combined -fiUyi(lS-l-2S-|-3S) in the SBS 
relative to the WBS is thus mainly due the y(lS). This 
supports the assertion made in the original work [7[ that 
the suppression level of the ground-state T is the most 
sensitive observable for color screening. In the SBS, re- 
generation contributions are small but rather significant 
especially for the Xb, helped by its large inelastic reaction 
rate. They decrease by a factor of 2 upon decreasing the 
bottom cross section from 3.2 /ib (PHENIX central value) 
to 1.6 /xb (STAR central value), since bottom states are 
in the canonical limit at RHIC. 



B. LHC Results 

Building upon previous Y calculations for the LHC in 
Pb-Pb collisions at y / s]v]v=5.5 TcV 0] we examine the 
ccntrality dependence at y / ]T/v"/v=2.76TcV, displayed in 
Fig. [8] for inclusive and direct T, as well as for direct 
T' and Xb- The higher temperatures at LHC (relative 
to RHIC) reinforce the stronger suppression in the weak- 
binding scenario compared to strong bindin g, f or inclu- 
sive and direct T production. The CMS data [llj are now 
significantly better described in the SBS, suggesting little 
effect of color-screening on the ground state even under 
LHC conditions. The initially (would-be) produced ex- 
cited states exhibit close to complete suppression in both 
scenarios, but significant regeneration occurs in the SBS, 
due to a much larger equilibrium limit. This should be a 
measurable effect, especially in connection with p t spec- 
tra (not discussed in our present work). In addition, the 
feeddown of the regenerated excited states accounts for a 
significant portion of the observed (inclusive) T(1S) re- 



generative component. Without feeddown, the regenera- 
tion effect is ~5%, while including feeddown it is roughly 
twice as large, ~10%. The understanding of the behavior 
of the higher T' and Xb states is thus essential to a proper 
interpretation of the observed T ground-state 1 production 
at LHC. 



C. Comparison to Other Recent Work 

In this section we briefly discuss the results from other 
recent calculations of Y production in URHICs in com- 
parison to our findings. 

In Ref. @ Y suppression at RHIC was studied using 
gluo-dissociation with vacuum binding energies, folded 
over an ideal hydrodynamic evolution. Dissociation 
temperatures were implemented according to potential- 
model results based on V or U potentials. In line with 
our results for the SBS, the y(lS) remains essentially 
unsuppressed in Au-Au(y / s]v r /v=0.2 TeV) collisions, while 
the yields and p^-spectra of excited states show signifi- 
cant sensitivity to variations in the dissociation tempera- 
ture. The larger dissociation temperatures following from 
V=U give results close to our SBS. 

In Ref. [lH Y suppression was computed using an in- 
medium Cornell potential (l8| including a perturbative 
imaginary part and the effect of anisotropies. The imag- 
inary part corresponds to quasifree dissociation (22l . |23| , 
and was folded over an "extended anisotropic hydrody- 
namics" evolution. The resulting suppression in Pb-Pb 
at LHC is somewhat weaker than in our WBS, partly be- 
cause the anisotropies decrease the Y dissociation rates. 

In Ref. [l4| the gluo-dissociation cross section was 
extended to approximately incorporate the string term 
in the bb potential. With in-medium binding energies 
but with a rather small coupling constant of a s =0.2, a 
schematic estimate of the T suppression in central Pb-Pb 
at LHC is roughly compatible with CMS data. 

Close to completion of our work, Ref. [45| appeared 
where bottomonium production at RHIC and LHC is 
computed within the same rate-equation approach 0, H(| 
as employed here, including regeneration. The Y dissoci- 
ation widths were taken from a NLO perturbative QCD 
calculation folded over a schematic viscous fireball evolu- 
tion with an ideal-gas equation of state. No constraints 
from lattice correlators on the Y spectral functions were 
evaluated. Contrary to our results, a preference for an in- 
medium binding scenario has been inferred and no signif- 
icant regeneration contributions arc found even at LHC. 

V. CONCLUSIONS 

We have studied bottomonium production in ultrarel- 
ativistic heavy-ion collisions at RHIC and LHC utiliz- 
ing a kinetic rate-equation approach in a quark-gluon 
plasma background, including suppression and regenera- 
tion mechanisms dictated by detailed balance. Our cal- 
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culations arc based on previous work in Ref. Q where, 
in particular, two scenarios for bottomonia in the QGP 
were evaluated, i.e. vacuum and in-medium bound-state 
properties (vis-a-vis strong and weak-binding scenarios). 
At the time, the latter seemed to be the theoretically bet- 
ter motivated case. The present work, however, suggests 
a re-evaluation of this view. Firstly, the strong-binding 
limit provides a more natural explanation of the near 
temperature independence of bottomonium correlators 
now available from thermal lattice QCD. Furthermore, 
the application to recent STAR and CMS data at RHIC 
and LHC (including updated open and hidden-bottom in- 
put cross sections) indicates that a more stable T ground 
state gives a better description of its nuclear modification 
factor. These findings corroborate the picture that the T 
ground state is not much affected by color screening up to 
fairly high temperatures of ~3-4T c . We also found that 
regeneration contributions, which have been neglected in 
most previous works, are not negligible, especially for 
strong binding at the LHC. 

Several lines of future studies emerge. The background 
medium needs to be improved using realistic hydrody- 



namic evolutions. Some work in this direction has al- 
ready been conducted. The regeneration term should be 
calculated with time-dependent bottom-quark spectra, 
constrained by B-meson observables. Based on those, 
transverse-momentum spectra and elliptic flow of bot- 
tomonia need to be calculated (again, pertinent work has 
begun). In all these, the consistency to the charmonium 
and open-charm sector should be maintained. A com- 
prehensive description of heavy-flavor observables, based 
on theory constrained by lattice QCD, remains the long- 
term goal, which appears to be achievable. 
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FIG. 7: The nuclear modification factor for T(1S+2S+3S) including e + e~ branching (top row, compared to STAR data [T3]), 
direct T (second row), T' (third row) and Xb (bottom row), as a function of centrality in Au-Au(y / ijvjv=0.2 TeV) collisions at 
RHIC. The left column corresponds to the weak-binding scenario, the right one to the strong-binding scenario. In each panel, 
CNM effects alone are shown by the green band, CNM plus QGP suppression by the blue band, regeneration by the dashed 
pink line and the total by the red band. If regeneration is small, red and blue bands largely overlap. 
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FIG. 8: The nuclear modification factor for inclusive T (top row, compared to CMS data [Tl[), direct T (second row), T' 
(third row) and \b (bottom row), as a function of centrality in Pb-Pb(y / ijviv=2.76 TeV) collisions at LHC. The left column 
corresponds to the weak-binding scenario, the right one to the strong-binding scenario. In each panel, CNM effects alone are 
shown by the green band, CNM plus QGP suppression by the blue band, regeneration by the dashed pink line and the total 
by the red band. 



